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ABSTRACT 

We present a stacking analysis of the complete sample of Early Type Galaxies (ETGs) in the Chandra 
COSMOS (C-COSMOS) survey, to explore the nature of the X-ray luminosity in the redshift and stellar 
luminosity ranges 0 < z < 1.5 and 10® < L^ILq < 10*^. Using established scaling relations, we subtract 
the contribution of X-ray binary populations, to estimate the combined emission of hot ISM and AGN. 

To discriminate between the relative importance of these two components, we (1) compare our results 
with the relation observed in the local universe Lx.gas for hot gaseous halos emission in ETGs, and 

(2) evaluate the spectral signature of each stacked bin. We find two regimes where the non-stellar X-ray 
emission is hard, consistent with AGN emission. Eirst, there is evidence of hard, absorbed X-ray emission 
in stacked bins including relatively high z (~ 1.2) ETGs with average high X-ray luminosity (Lx-lmxb ^ 

6 X 10"*^ erg/s). These luminosities are consistent with the presence of highly absorbed “hidden” AGNs 
in these ETGs, which are not visible in their optical-IR spectra and spectral energy distributions. Second, 
confirming the early indication from our C-COSMOS study of X-ray detected ETGs, we find significantly 
enhanced X-ray luminoaity in lower stellar mass ETGs (Lx < 10**Lq), relative to the local Lx^gas ‘x 
relation. The stacked spectra of these ETGs also suggest X-ray emission harder than expected from 
gaseous hot halos. This emission is consistent with inefficient accretion 10^^ - 10^'^MEdd onto Mbh ~ 
10®- 10*^Mo. 


Subject headings: galaxies: elliptical - surveys - X-rays: galaxies 


1. Introduction 

Diffuse X-ray emission from normal early type galaxies (ETG) was first detected in the Virgo cluster with the 
Einstein X-ray Observatory (Eorman et al. 1979). Although in a few cases (e.g., M86) the displacement of the X-ray 
emission from the stellar body of the galaxy pointed to hot gaseous halos, the ~ arcminute resolution of the Einstein 
IPC data could not disentangle observationally hot halos from the populations of low-mass X-ray binaries (LMXBs), 
which was likely to dominate the X-ray emission in a large number of the observed ETGs (Trinchieri & Eabbiano 
1985). Nonetheless, several authors discussed how the X-ray emission and the inferred hot gas contribution could be 
used to constrain both the binding mass of the galaxies and the physical evolution of the hot halos, including feedback 
from SNIa and AGNs, and interaction with intracluster and group media (see e.g., Eabbiano 1989; Canizares et al. 
1987; Ciotti et al. 1991; David et al. 1991; White & Sarazin 1991; Mathews & Brighenti 2003; Ciotti & Ostriker 
2007). 
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Thanks to the high spatial resolution of Chandra and imaging/spectral capabilities of the ACIS detector (Garmire 
et al. 2003), it has become possible to directly detect and characterize the different contributions to the X-ray emission 
of ETGs in the near universe. With Chandra, LMXB populations - both associated with Globular Clusters (GC) and 
in the stellar field of the ETG - have been studied to distances of several tens of megaparsecs (e.g., Sarazin et al. 2000; 
Kundu et al. 2002; see review by Eabbiano 2006). The integrated LMXB luminosity correlates well with the integrated 
ETG stellar luminosity (e.g., Kim & Eabbiano 2004; Gilfanov 2004). The scatter in this relation is a factor of 2-3, and 
it depends on both GC content (Kim & Eabbiano 2004), and the age of the stellar population (Kim et al. 2009; Kim & 
Eabbiano 2010; Boroson et al. 2011; Lehmer et al. 2014). These detailed high resolution studies allowed Boroson et al. 
(2011, BKEl 1) to accurately estimate the gaseous contribution to the X-ray luminosity Lx,gas of a sample of 30 ETGs 
within 32 Mpc. After subtracting LMXB, stellar, and nuclear contributions to the X-ray luminosity, BKEl 1 found 
steep scaling relations between Lx,gas and both the stellar Lx luminosity and the hot gas temperature T . Subsequent 
work by Kim & Eabbiano (2013, KE13; 2015, KE15) suggests that for Lx,gas > 10^° erg s“' the hot gas around ETGs 
is virialized, giving a new way to measure the total mass of these galaxies. 

Beyond the local universe, Chandra Deep Eield studies (both detections and stacking of almost 1,000 ETGs) 
found that the mean Lx of Lb = lO'^L© ETGs remains equal or mildly increases over z = 0 - 0.7 (Lehmer et al. 2007; 
Tzanavaris & Georgantopoulos 2008) and that the ratio of the soft X-ray luminosity to the B-band luminosity evolves 
mildly as E (0 5-2 keV)/^^ _|.^y.i±o.7 ^ ~ p2 (Danielson et al. 2012), suggesting some heating mechanism 

preventing the hot gas from cooling, consistent with mechanical heating from radio AGNs. The evolution of the X-ray 
binary contribution to the X-ray emission (both of LMXBs and high mass binaries) across cosmic time, modeled by 
Eragos et al. (2013a), led Eragos et al. (2013b) to conclude that energy feedback from X-ray binaries is also important 
in the early universe. 

The COSMOS survey, covering ~ 2 square degrees of sky with multi-wavelength coverage, which includes 
Chandra X-ray coverage (Scoville et al. 2007; Elvis et al. 2009), has increased the available sample of ETGs to over 
6,600 galaxies Civano et al. (2014, C14). This survey provides a good z coverage out to z ~ 2.5 (and a few galaxies 
are detected out to z ~ 5). C14 studied the sample of 69 X-ray detected ETGs in the Chandra-COSMOS survey 
(C-COSMOS, Elvis et al. 2009), covering aredshift range 0 < z < 1.5, to explore redshift dependencies in the BKEll, 
KE13 scaling relations. In particular, C14 focused on the relation between Lx,gas and the K-band stellar luminosity 
{Lx) of ETGs which is a good proxy of stellar mass for old stellar systems. The latter has been evaluated in the local 
Universe by BKEll as Lx,gas for sources with Lx,gas < 2 x 10^' erg s“', while including sources with 

Lx,gas up to 9 X 10'*^ erg s“' KE13 found a steeper Lx^gas ^ relation. C14 found that this latter relation also 

holds for X-ray detected ETGs with Lx,gas ^ xlO"^^ erg s ', suggesting that the hot gas in these more distant galaxies 
may be virialized as well. However, more distant (z > 0.9) and younger stellar age galaxies generally fall outside the 
relation observed in the local sample, and are X-ray overluminous with respect to their Lx- Eor these sources C14 
suggested that the higher X-ray luminosity may be due to the presence of hidden AGNs (as suggested by observed 
hardness ratios), merging phenomena, which are expected to increase the halo X-ray luminosity, or nuclear accretion 
(that could also induce AGN activity during the merging phase, Cox et al. 2006). 

In this work we perform a stacking analysis of the C-COSMOS ETG data in order to investigate the nature of 
fainter. X-ray undetected ETGs. We compare their properties with those of the sample of local ETGs investigated by 
BKEll and KE13, and with the sample of ETGs detected in the C-COSMOS survey studied by C14. This stacking 
study also allow us to further explore redshift dependencies. Our aim is to explore the properties of the non-stellar X- 
ray emission, i.e. the X-ray emission resulting from the subtraction of the expected X-ray binary contribution (taking 
into account the cosmic evolution of the X-ray binary population Eragos et al. 2013b). The resulting emission will 
probe both hot X-ray halos and the possible emission of active nuclei in these galaxies. In particular, the large number 
of low-stellar-mass ETGs in the sample will provide a unique probe of the existence of relatively low-mass massive 
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black holes in their nuclei. These black holes should produce low-luminosity nuclear X-ray sources, if they are fueled 
by the ETG stellar outgassing (Volonteri et al. 2011). 

The paper is organized as follows: Section 2 describes the sample selection and the data sets used in this work, 
Section 3.3 is dedicated to the description of the adopted stacking procedure, in Section 4 we discuss our results, and 
Section 5 is dedicated to our conclusions.In the following, we adopt the standard flat cosmology with Qa = 0.73 and 
Hq - 70 km s“' Mpc“*. 


2. Sample Selection 

The COSMOS survey (Scoville et al. 2007) provides a multi-wavelength characterization of over 1.5 million 
galaxies at redshift extending out to 5 in 2 deg^ of the sky. The Chandra X-ray coverage (C-COSMOS Elvis et al. 
2009) over 0.9 deg^, gives us the means to study X-ray emission of ETGs, extending local studies to higher redshifts, 
and explore the redshift evolution of hot halos with large samples of galaxies. 

Our ETG sample is selected from the most recent version of the COSMOS photometric catalog by Ilbert et al. 
(2009). The spectral energy distribution (SED) identification performed by Ilbert et al. (2009) makes use of 7 elliptical 
galaxy templates with ages from 2 to 13 Gyr (Polletta et al. 2007). The SED templates were generated with the 
stellar population synthesis package developed by Bruzual & Chariot (2003), assuming an initial mass function from 
Chabrier (2003) and an exponentially declining star formation history. We also collected additional galaxy properties 
(stellar mass, age and star formation rate) from Ilbert et al. (2010), and we made use of the most recent near-infrared 
photometry from the Ultra Deep Survey with the VISTA telescope (Ultra-VISTA McCracken et al. 2012)'. 

We selected galaxies classified as E or SO and that have a reliable photometric redshift (Ilbert et al. 2009). In 
addition we selected only sources with Ultra-VISTA detection both in J and K bands, since these are needed to evaluate 
the rest frame Lk (see below). This narrows our sample to 10,808 (the COSMOS ETG sample considered in C14). 
We then selected only the 6436 ETGs falling in the C-COSMOS field, and then excluded the 48 X-ray detected ETGs 
in common with the C-COSMOS identification catalog (Civano et al. 2012). This leads us to our final sample of 6388 
C-COSMOS ETGs - not detected in X-rays - which constitutes the starting point of the stacking analysis described in 
this paper. Note that we selected sources with photometric redshift estimate in the catalog by Ilbert et al. (2009), while 
C14 selected their sources to have a spectroscopic or photometric redshift estimates in the C-COSMOS identification 
catalog, so the sample of ETGs selected by the latter authors (~ 6600 sources) is not fully contained in our sample. 

Eollowing C14 we evaluated the K band luminosity of our sources from the Ultra-VISTA K-band aperture mag¬ 
nitude of the COSMOS photometric catalog. To calculate the rest frame K band luminosities we assumed a power 
law spectrum (as appropriate over this small frequency range) /y cx v'® with a - { J - K)I2.5 log (vjIvk) (where J and 
K are the AB magnitudes from the COSMOS photometric catalog and vj and vk are the corresponding frequencies), 
and express the rest frame luminosities in solar luminosity as x (1 -i- z)"^' x (D^/10)^, where z 

is the photometric redshift from the COSMOS catalog, Di is the luminosity distance expressed in parsecs, and ext is 
the aperture correction parameter computed from the EWHM measured for all the COSMOS extended sources in the 
Hubble ACS data. 

We compared the distributions of the physical parameters of the 6388 C-COSMOS ETGs with the entire 10,808 


'Note that stellar mass reported in COSMOS photometric catalog, and obtained from the Bruzual & Chariot (2003) best-fit templates, do not 
incorporate the Ultra-VISTA K-band NIR data. Although this is only an issue if there is any residual star formation in the ETGs, Lk and M* may 
not necessarily be completely equivalent. 



-4- 


COSMOS ETG sample finding no significant difference, demonstrating as expected that the sample of ETGs not 
detected in X-ray is representative of COSMOS ETGs overall. Note that the luminosity threshold is a function of 
redshift because of the uniform magnitude limit of the Ultra-VISTA COSMOS survey (K^ ~ 24, McCracken et al. 

2012) . The specific star formation rate for the majority of the ETGs is in the regime where Ilbert et al. (2010) defines 
galaxies as quiescent, although higher star formation rate is observed in young ETGs. The properties (luminosity, 
mass, star formation rate and age) of the COSMOS ETGs are in agreement with those reported by (Moresco et al. 

2013) for a sample of ETGs selected using multiple criteria (including color-color diagrams, spectra, star formation 
rate and SED classification). Masses and luminosities are in agreement with those of the K-selected elliptical galaxies 
of Jones et al. (2014), extracted from a small area of the COSMOS field. 


3. Analysis and Results 

Below we describe the procedures followed in our stacking analysis. The results are summarized in Tables 1 and 

2 . 


3.1. Count Extraction at ETG Positions 

We used all 49 C-COSMOS observations to evaluate source and background counts for our sample of ETGs, 
excluding data in the ACIS-I chip gaps. We also excluded data from X-ray detected sources removing counts within 
a 10” radius of all sources in the C-COSMOS source identification catalog (Civano et al. 2012). We replaced the 
excluded counts by sampling the Poisson distribution of the pixel values nearby background regions with dmfilth 
CIAO task^. When a source lay in the field of view of more than one observation, we combined the data from the 
different observations using the prescriptions from Puccetti et al. (2009), excluding observations where the source was 
more than 10 acrmin off-axis from the aim point. 

Due to the variation of the Chandra Point Spread Eunction (PSE) radius Rpsp^, we used an extraction radius 
that was determined by the source position in each observation in order to extract an optimal fraction of X-ray photons 
from the diffuse halo. Eor sources with morphological information available in the HST -ACS catalog (Leauthaud et 
al. 2007) we adopted the source semi-major axes R as a measure of its physical extent. Eor sources with R > Rpsp 

we adopted values of that was the convolution of R and Rpsp -^jiR^ + otherwise we adopted re^t = Rpsf 

as we did for sources without morphological information in HST-ACS catalog. We notice however that sources with 
R > RpsF represent only ~ 1% of our sample, with extraction radii typically ranging between 1” and 7”. The average 
exposure time for each source within the extraction circle, E, was then evaluated. The histogram of extracted ETG 
position counts is shown in Eigure 1 (blue line). The small fraction of extraction regions with counts above ~ 5 
represent off-axis sources with large extraction radii. 


^http://cxc.harvard.edu/ciao/ahelp/dmfi1th.html 

^In the following we consider the Chandra-AClS radius enclosing 90% of the energy. 




-5- 


3.2. Background Evaluation 

To evaluate the background we extracted counts from an annulus centered on the source coordinates, with inner- 
radius 2r„, and outer radius 4r„,. Due to the very low signal, the background represents the largest source of uncer¬ 
tainty in this kind of analysis. To estimate the systematic uncertainties introduced by the evaluation of the background, 
we compared the results of three different methods, obtaining consistent results. First we simply evaluated the back¬ 
ground contribution to the total source counts from the mean count value of the pixels in the whole background 
annulus. Second, we adopted the clipping procedure of Willott (2011, method b), which evaluates the mean count 
value fj - and the standard deviation cr - of all pixels within a circle of radius 4r„, around each source position and then 
sets all pixels with count value larger than ju + 3 (t to ji. This method suppresses the noise in the background annulus, 
but also removes signal from the source"^. Lastly, we divided the background annulus into eight 45 degree sectors, 
and excluded all the sectors in which the mean count value was larger than jU -i- 3cr, thus excluding regions of higher 
background without removing signal from the source. The net counts N obtained with all three methods are consistent 
within statistics, so in the following we adopted the clipping procedure of Willott (2011, method b). 

To test the robustness of our results we also extracted counts from positions on the sky offset from the ETG 
coordinates by random distances between 15” and 30” both in right ascension and declination. The net counts of the 
non-source positions are shown in Figure 1. They follow an approximatively Gaussian distribution centered on zero, as 
expected from background fluctuations, while the ETG counts show an excess of positive signal. To test if our results 
were dominated by a sub-population of ETGs we repeated the stacking procedure on sub-samples randomly drawn 
from the complete sample, progressively halving the subsample, obtaining results compatible with the full sample. 


3.3. Stacking Procedure 

To evaluate the count rate C for each stacked bin we stacked the observed source net counts in 0.5-7 keV range 
N and exposure E to obtain C - 2,W,7 2, £/. Erom these counts we subtracted the expected contributions from 
stellar sources (LMXBs and Ultra Luminous X-ray sources, ULXs), and evaluated the average bin rest-frame X-ray 
luminosity as explained below. We corrected Lk to take into account the evolution of the stellar population (i.e., the 
fading effect, see C14). Note that, at a given redshift, younger galaxies have larger fading correction than older ones. 
To compare our results with those of KE13 and C14 we considered sources with 9 < log (Lk) < 13. 

We constructed the stacks by binning in both Lk and z, to ensure that the binning scheme did not produce any 
significant bias on the results. Eor the Lk binning, we first divided our sample in four log (Lk) bins, 9 - 10, 10 - 11, 
11-12 and 12 - 13. We then ordered the sources in each bin in increasing redshift, and progressively stacked them 
until we reached a minimum signal to noise ratio of 3; if this 3cr threshold could not be reached we stacked the signal 
from all the remaining sources in the bin. Examples of 0.5-7 keV images for these stacking bins are presented in 
Eigure 2. Eor the z binning, we followed the same procedure, dividing the sample into three bins 0 - 0.5, 0.5 - 1 and 
1 - 1.5, and then ordered the sources in each bin in increasing Lk- 


'^For the rare cases where ji + 3cT < 1, only pixels with two or more counts were clipped. 
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3.4. Evaluation of Non-stellar Lx for Each Stacked Bin 

To subtract the LMXB contribution from the total X-ray emission, we evaluated the average LMXB contribution 
to each bin using the 0.3-8 keV relation from Fragos et al. (2013b), as in C14; 

log(Lz.MVB/M.) = 30.259 - 1.505 x log(age) - 0.421 x log^ (age) -H 0.425 x log^ (age) -H 0.135 x log'* (age), (1) 

where Lx is expressed in erg s“', M, is the average stellar mass of the stacking bin expressed in solar masses and age 
is the average stellar age of the stacking bin expressed in Gyr. ^ We also evaluated the contribution to X-ray luminosity 
from ULXs that could be present in the galaxies using the scaling relations from Calzetti et al. (2007) and Mineo et 
al. (2012) between the ULX luminosity and the galaxy star formation rate. ULX contribution is more than an order of 
magnitude smaller than Lx-lmxb even for the least luminous Lx bins. 

We converted the average 0.3 - 8 keV LMXB luminosities in each bin to count rates in the 0.5 - 7 keV interval 
using a power law model with index F = 1.8 (consistent with the typical spectrum of LMXBs) and Galactic absorption 
Nh - 2.6 X 10^° cm“^. For these parameters the conversion factor used from Chandra count rates to fluxes is 1.04 x 
10^** erg cm“^ s“' cts^*®. We then subtracted these count rates from the stacked signals and converted back the 
resulting count rates into luminosities using a model comprising Galactic Nh plus a thermal APEC component with 
temperature depending on the source total X-ray luminosity. We adopted kT - 0.7, 1 and 2 keV for X-ray luminosities 
of Lx < 10“**, 10"** < Lx < 10"*^ and Lx > 10*^ erg s“', respectively, following BKFll and Dai et al. (2007). Finally, 
we K-corrected the luminosity to the 0.3 - 8 keV rest frame band. These LMXB-subtracted luminosities, that we 
indicate as Lx-lmxb^ ure expected to include the emission from the hot gaseous halos, but may also include a possible 
contribution from any AGN harbored within these ETGs. 

Eigure 3 shows the results in the Lx-lmxb - Lx plane for both Lx and z binning (Left and Right panels, respec¬ 
tively), where on the upper axes we plot M* as taken from COSMOS catalog. In both binning schemes the overall 
scatter plots are consistent, showing that we do not suffer from strong binning biases. By construction, the L/^-first 
binning provides information on the extreme redshift bins (both high and low), while the z-first binning covers more 
fully the Lx range. The bins are numbered in Eigure 3, following Tables 1 and 2. We use different colors to visually 
divide the stacking bins into four groups: red, for bins X-ray over-luminous with respect to the local Lx-lmxb and Lx 
relation (BEKll, KED, see C14); yellow, for X-ray under-luminous bins; and grey and green for bins following the 
local relation with Lx-lmxb lower and higher than 10"** erg s *, respectively. 

As noted in Section 1, Eq. 1 has an observed spread of a factor 2-3, which can be related to both the presence 
of GC LMXB in the ETGs and to the spread in age of the stellar population (Kim & Eabbiano 2004, 2010; BKEll; 
Lehmer et al. 2014). Given the parameter space of Eigure 3, this spread will not change our conclusions. 


3.5. Hardness Ratios 

Hardness ratios (HRs) were evaluated to characterize the ETG average spectral properties. HR - (H-S)f(H+S), 
where S is the soft 0.5 - 2 keV counts and H the hard 2-7 keV counts. These HRs are those of the LMXB-subtracted 
emission. The contributions of these sources were evaluated and subtracted by applying the procedure in Section 3.4 
to each energy band. We evaluated the errors on HRs as due to both uncertainties on counts and to uncertainties on 
estimate of LMXB contribution due to spread of stellar age and mass in the bins. 


^These average quantities, as well as their standard deviation in each bin is reported in Tables 1 and 2. 
®http: //cxc. harvard. edu/toolkit/pimnis. j sp 
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In Figure 4 we show the HRs for the stacking bins using the same color code as in Figure 3. We plot HR versus 
z in the Left panel, using the L/^-first binning, which gives a better z sampling. We plot HR versus Lk in the Right 
panel, using the z-first binning, to explore the Lx behavior. The left panel of Figure 4 shows a distinct clustering in 
HR-z space of bins from separate loci of the Lx-lmxb - Lx plane, as shown by the clusters of different color points. 
The Right panel shows that ETGs with larger stellar mass (Lx) tend to have lower HR. Given the scaling relations of 
ETGs (BKEll; KE13; KE15), this could just reflect a relation with total binding mass, with more massive galaxies 
being more efficient at retaining large hot gaseous halos. 


4. Discussion 

Eigure 5 summarizes our results in the Lx-lmxb - Lx plane and compares them with the COSMOS ETG X-ray 
detections (C14) and with the near universe ETG data (BKEl 1, KE13). We show the local KE13 relation with a solid 
line, and the scatter in this relation with the two dashed lines. These lines are chosen to pass through the points for 
the most Lx^gas luminous source (M87) and the most Lx luminous source (NGC1316), respectively. The COSMOS 
stacked-ETG results are overall consistent with the C14 COSMOS detections, but significantly increase the coverage 
at lower total stellar luminosities / masses. While most of the detected signals are compatible with the local KE13 
relation, there are several detections that lie both to the left and to the right side of it, that is, bins that appear to be 
respectively over and under-luminous in X-ray with respect to their Lx or M». 

In Eigure 5 we also show the hardness ratios for each stacked bin with lines with slopes proportional to HR. We 
find positive HRs (i.e., harder spectra) in stacked bins outside the dashed line on the left. These are the lower Lx 
galaxies with significant X-ray excess (the red points in Eigure 4, see also C14). We also find hard HR > 0 in some 
points with Lx-lmxb > 10^* erg s ' and z > 1 (the three top green points in Eigure 4), to the left of the solid line. 

Instead, lower HR ratios (i.e., softer spectra) are found in bins following the local relation (between the two 
dashed lines). This is especially so for those to the right of the solid line. This suggests that they have a larger than 
average hot gas content up to z ~ 1. 

Eigure 6 compares the HR-z scatter diagram with a set of spectral models. On the Left, we show the expected 
HRs for thermal APEC (dashed lines; from bottom to top, kT — 0.7, 1, 2 and 3 keV) and power-law models (solid 
lines; fixed slope F = 2 and increasing intrinsic absorption from bottom to top, Nh - 0, 10 , 10 and 10 cm ). 
The APEC models are the expected emission of thermal gaseous emission, while the power-law models are typical 
AGN spectra. To illustrate cases of mixed emission from a hot gas component and a heavily obscured AGN, on the 
right panel of Figure 6 we plot composite models consisting of kT =1 keV APEC emission with line of sight Nh plus 
power-law models with slope F = 2 and intrinsic absorption Nh - 10^^ cm“^. From the bottom to the top the relative 
normalization of the two models varies, to reflect an increasing contribution of the power-law model to the total flux, 
ranging from ~ 0.2% to ~ 100%. 

The X-ray over-luminous bins (red points) have generally high HRs consistent with being dominated by AGN 
emission, either with “normal” AGNs with absorbing columns > 10^^ cm“^, or in the case of the composite gas -i- 
AGN emission even highly obscured Compton Thick AGNs {Nh > lO^'' cm“^, Risaliti et al. 1999; Levenson et al. 
2006; Lehmer et al. 2010). 

The bins following the local relation with Lx-lmxb < 10“" erg s * (grey points) show HR ~ (-0.5) - 0.3, lower 
than those of the over-luminous X-ray bins. A soft HR first suggests emission from hot gas. However, Figure 6 shows 
that pure APEC thermal models with solar metallicity always have HR < -0.6. Hence these objects have HR values 
inconsistent with pure gaseous emission. 
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These HRs can instead be explained with combinations of thermal and AGN emission. Figure 6 Right shows 
models with varied ratios of APEC and power-law emission. In this model the observed HRs imply AGN components 
accounting for 15% - 60% of the X-ray luminosity. Given the X-ray luminosities of this sub-sample (Figures 3, 5), 
these ratios then imply AGN luminosities ranging from ~ 5 x 10^® erg s“' to ~ 5 x erg s“'. This range is in 
agreement with the low-luminosity Chandra X-ray sources in the near universe sample (BKFll, Pellegrini 2005), 
which have luminosities ranging from ~ 10^^ erg s“' to ~ 5 x 10"^* erg s“'. 

The bins with Lx-lmxb > lO'** erg s ' (green points) have hard HR > 0 incompatible with a low absorption 
power-law model (e.g., bins 21, 22 and 25). The three highest redshift bins in this group (e.g., bins 28, 29 and 30) 
have even harder HR > 0.5 that are compatible only with absorbing columns larger than 10^^ cm“^. This suggests that 
these galaxies harbor (on average) highly obscured Compton Thick AGNs that dominate their X-ray luminosity. 

In summary, our stacking analysis suggests that on average AGN emission is likely to be present in ETGs. These 
AGNs may be faint and hidden in hot gas dominated ETGs, which are consistent with the near Universe scaling 
relation of virialized halos (KF13, KF15, see C14). In contrast, AGN emission appears to dominate the non-stellar 
X-ray emission in two classes of ETGs: (1) low-stellar mass ETGs (log (LkILq) < 10.5), and (2) higher redshift ETGs 
(z > 0.5, >1) especially those overluminous in X-rays relative to the near universe Lx-lmxb - Lx scaling relation. 
The latter group has HR values suggesting highly absorbed Compton Thick AGNs. We note that these AGNs were not 
visible in the COSMOS colors and spectra, since AGNs were excluded in our ETG sample selection process (C14). 


4.1. Accretion Regimes for the AGNs 

Figure 7 shows the results of the L/f-first stacking, with Lx-lmxb plotted against stellar mass M, and the mass of 
the nuclear black hole, Mbh, obtained from the Mbh - M, relation of Kormendy & Ho (2013). Evaluating Mbh from 
the Mbh - M, relation McConnell & Ma (2013) would yield black hole masses between 0.1 and 0.3 dex smaller than 
the former. Similarly, one could derive Mbh from the Mbh - Lx relation from Graham (2007), obtaining black hole 
masses ~ 0.5 dex larger than the former. Uncertainties in these estimate are typically of a factor of ~ 2 - 3 (Tremaine 
et al. 2002; Graham 2007) and evolution to z ~ 1-2 may provide a systematic black hole mass increase of a similar 
factor (Merloni et al. 2010). These effects will not change our conclusions in the 3-dex wide range covered by Figure 
7. 


On the right hand axes we plot Lhoi - IQ Lx-lmxb, an appropriate bolometric correction for low-luminosity AGNs 
that radiate at a few 10s % of Lxdd (Kollmeier et al. 2006; Steinhardt & Elvis 2010). Given this nominal Lhoi we can 
plot diagonal lines of constant LhoilLEdd- These lines range from 10“^ to 10“® LtoilLEdd- 

We see that the massive high Lx-lmxb galaxies (green points) and the low-mass X-ray overluminous (red points) 
galaxies have implied Eddington ratios in the 10^^ to 10^^ range. RIAFs are generally believed to apply to all black 
hole accretion below IQ^^LjLEdd, so we would expect all these sources to be RIAFs. If the bolometric correction for 
RIAFs is smaller than that considered here, as may well be the case (Narayan & Yi 1995), the Lhoi values on the y-axis 
should be correspondingly lower (as the Eddington ratios on the currently plotted dashed lines). 

Alternatively the true Lhoi of the three high z bins of massive high Lx-lmxb galaxies (green points 28, 29, 30) 
could be greatly underestimated. The HRs of these bins also imply highly absorbed emission {Nh > 10^^ cm“^. 
Figure 6). If these are then Compton Thick AGN they may have much larger intrinsic luminosities, with factors of 
100 being quite plausible (e.g., Risaliti et al. 1999; Levenson et al. 2006). They would have intrinsic luminosities 
A) 5-10 keV ~ ^ LhoilLEdd in the ~ 10^^ range, so that they could lie within the normal AGN regime. 

To investigate this possibility, we built SEDs of ETGs in our sample using optical-infrared magnitudes reported 
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in COSMOS photometric catalog, since a “hidden” AON may give rise to an IR excess with respect to an average ETG 
SED. We then fitted these SEDs with the same elliptical galaxy template used by Ilbert et al. (2009) without finding 
any significant deviation. To confirm this, we evaluated the E 2 - 10 -LMXB luminosity of the stacking bins assuming a 
power-law spectrum (as appropriate for AON spectra) and compared that with their IR luminosity. In all cases the X- 
ray luminosities are at least 3 order of magnitudes lower than the IR luminosities, even assuming absorbing columns 
~ 10^^ cm“^, indicating that at IR wavelength a putative AON component is expected to be significantly smaller than 
the galactic emission. 

The galaxies that are consistent with the local Lk - Lx-lmxb gaseous relation (grey points), with X-ray colors 
suggesting thermal emission plus a possible AGN component, have smaller Eddington ratios ~ 10^^ to 10 ® (these 
ratios of course assume that the entire emission is due to an AGN; as shown in Eigure 6 , the AGN component may 
range from ~ 15 - 60%). 

We also plot in Eigure 7 the predictions from Volonteri et al. (2011) of Lx^acn for three values of M* for both 
efficient and radiatively inefficient accretion flows (RIAEs Narayan & Yi 1995). These authors’ models use the stellar 
density profile scaling with stellar mass, and predict the accretion rate from the normal stellar mass loss due to stellar 
evolution, and the consequent massive black hole luminosity for various radiative efficiencies^. The X-ray luminosity 
coming from accretion is computed as Lx^gn!^E dd - ? 7 x? 7 'h, where tjx - 0-I is the X-ray fraction of the bolometric 
luminosity, and m - OAMjMzdd is the accretion rate in Eddington units with a 10% accretion efficiency. Radiatively 
efficient accretion flows have rj - while inefficient ones have 77 = mlihcr, and for them m < ihcr = 3 x 10“^. 

In order to compare our results with Volonteri et al. (201 1) models we converted the 0.5 -7 keV LMXB subtracted 
counts into 0.5 - 10 keV band luminosities, ^ jg assuming a pow-law spectrum as appropriate for AGN spectra. 
Based on these models, inefficient accretion into the nuclei could explain even all the Lx-lmxb of the massive, high 
Lx-lmxb ETGs (green points), except for the most luminous ones, which are likely to be Compton thick AGNs (see 
above). The low-mass X-ray overluminous ETGs (red points) are also consistent with the Volonteri et al. model. 
Considering that a large fraction of the X-ray emission is due to hot gas (see above), the grey points show instead a 
possible AGN emission even lower than predicted for radiatively inefficient accretion. In these ETGs, the Lx-lmxb is 
consistent with the local Lx,gas - Lx relation. 

This result can have the following possible explanations, that are not mutually exclusive. 

1. The Volonteri et al.’s models assume that the whole of their estimated mass accretion rate (coming from the 
instantaneous mass loss of all stars, and obtained by integrating over the whole stellar density profile) indeed 
reaches the massive black hole, while this represents an upper limit on the mass accretion rate of the massive 
black hole (some gas may be outflowing, or may reach the massive black hole with a delay, or may possess 
angular momentum and form a disk, etc.). 

2. The accreting mass, that has already entered the accretion radius, is reduced on its way to the massive black 
hole, due to a wind from a RIAE (Narayan & Yi 1995; Blandford & Begelman 1999); similarly, close to the 
nucleus, there may be a heating mechanism for the accreting gas, as a nuclear jet, again impeding or lowering 
the accretion rate (e.g. Pellegrini et al. 2012; Yuan & Li 201 1). 

3. Accretion onto the black hole is (temporarily) not taking place, as may be the case during feedback-regulated 
activity cycles. Numerical simulations show that, during these cycles, AGN outbursts are followed by major 
degassing of the circumnuclear region, with a precipitous drop of the nuclear accretion rate; the outbursts are 


^We note that in Volonteri et al. (2011) models assume for ETGs a stellar age of 12 Gyr. We corrected for the average age 10® ^ Gyr of 
C-COSMOS ETGs, which yields slightly smaller emitted X-ray luminosities (see Eq. 9 in Volonteri et al. 2011 ). 
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separated by long time intervals during which the galaxy is replenished again by gas from the stellar mass losses, 
until a new nuclear outburst takes place (Ciotti et al. 2010). 

4. rjx might be lower than 1 as assumed here (see also Volonteri et al. 2011), or that tj at low m decreases faster 
than adopted in the Volonteri et al.’s models, so that lower Lx^gn should be predicted (e.g.. Yuan & Narayan 
2014). 

Our results extend to higher X-ray luminosities a recent analysis of ETGs in the AMUSE survey of Virgo and 
field galaxies with Chandra (Miller et al. 2015). AMUSE results indicate a significant presence of supermassive black 
holes in low stellar mass galaxies. Due to the lower redshifts of the AMUSE ETGs, that survey is sensitive to X-ray 
luminosities (~ 1.3x10^^ erg/s) one order of magnitude fainter are than our faintest bin 1 (~ 2.4x10^^ erg/s). Eor these 
dim sources the stellar wind accretion proposed by Volonteri et al. (2011) would predict higher X-ray luminosities than 
those observed. Miller et al. conclude that these black holes are probably powered by inefficient advection dominated 
accretion flows or, alternatively, by an outflow/jet component. We found that low-luminosity. X-ray “overluminous” 
ETGs can be adequately described by the RIAE fueled by stellar winds as proposed by Volonteri et al. (201 1). Instead, 
ETGs with Lx^lmxb > 10"^^ erg s“^ may host Compton Thick AGNs. Similarly, when comparing with the results of 
Lehmer et al. (2007) and Danielson et al. (2012) on Chandra Deep Fields, we extend to higher X-ray luminosities 
and to lower optical luminosities (that is, lower stellar masses), finding sources that significantly exceed the expected 
emission from hot-gas dominated ETGs, with log(Lj,g j !Lq) > 30 (see Figure 7 in Danielson et al. 

2012 ). 

We also note that, as reported by Comastri et al. (2015), the revision of black hole scaling relations by Kormendy 
& Ho (2013) indicate that the local mass density in black holes should be five times higher than previous estimates. 
These authors advance the possibility that a sizeable population of highly obscured AGNs with covering fraction of 
obscuring material ~ 47r would explain the revised upward value of the SMBH local mass density without exceeding 
the limits imposed by the X-ray and IR background. An example of these obscured AGNs can be that observed in the 
ETG ESO565-G019 with Suzaku and 5wi/f/BAT (Gandhi et al. 2013). 

Finally, we note that the hard X-ray slopes of the more luminous ETG bins could be due to contamination from 
X-ray jets. However, the combined 1.4GHz VLA COSMOS Large and Deep Project catalogs (Schinnerer et al. 2010) 
offer radio counterparts only to ~ 2% of our sample of ETGs. The presence of radio emitters among non X-ray 
detected ETGs will be addressed after the release of the JVLA-COSMOS Project data. 


5. Summary and Conclusions 

We have followed up our previous work (C14) investigating the X-ray / K-band properties of the 69 ETGs detected 
in the Chandra COSMOS survey (Elvis et al. 2009), with a luminosity and hardness ratio stacking analysis of the 
Chandra data for the sample of 6388 early type galaxies (ETGs) not detected in X-rays. The stacking analysis allows 
the investigation of the average features of individually X-ray undetected ETGs. Our ETG sample was selected as 
representative of normal galaxies, with no evidence of AGN emission in their spectra and multi-band photometry 
(C14, Section 2). Our purpose was to investigate the properties of the X-ray emission from hot gas and from possible 
“hidden” AGN components that may be revealed by their X-ray emission. To this purpose, the expected contribution 
of LMXBs was subtracted from the X-ray luminosity of each stacked bin, using established scaling relations with the 
K-band luminosity (BKFl 1). All luminosities were corrected to the local rest frame; X-ray hardness ratios were based 
on the observed counts (after subtraction of the LMXB contribution), and were compared to emission models for the 
appropriate average bin redshift. 


This analysis has led to the following conclusions; 
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1. On average AGN emission may be present in all ETGs. 

2. AGNs may be faint and hidden in hot gas dominated ETGs, which are consistent with the local universe scaling 
relation of virialized halos (KE13, KE15, see C14). 

3. AGN emission appears to dominate the non-stellar X-ray emission in low-stellar mass ETGs (log(L/i-) < 10.5 
in solar units), which show marked X-ray luminosity excesses relative to their K-band luminosity (a proxy of 
stellar mass), when compared with the local ETG Lx^cus^k relation of the local universe (BKEl 1). This result 
confirms the conclusions of C14 and extends them to lower stellar mass galaxies. 

4. AGN emission is also prominent in higher redshift ETGs (z > 0.5, > 1) especially those with X-ray luminosity 
> 10"^^ erg s“^. The latter group have HR values suggesting highly absorbed Compton Thick AGNs. We note 
that these AGNs were not visible in the COSMOS colors and spectra, since AGNs were excluded by our ETG 
sample (C14). 

5. Given the nuclear BH mass predicted by the Mbh-Lk relation (Kormendy & Ho 2013), the hidden AGNs in the 
X-ray over-luminous low stellar mass ETGs are consistent with the presence of radiatively inefficient accretion 
fueled by stellar outgassing (Volonteri et al. 2011). The highly absorbed higher Lx-lmxb AGNs hidden in 
massive, higher z ETGs instead suggest efficient accretion in highly obscured nuclei. 
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Fig. 1.— Distribution of net counts for randomly distributed extraction regions (red line) and for the selected sample 
of sources (blue line). 
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Fig. 2.— Examples of 0.5-7 keV images of the Lk selected stacking bins. 
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Table 1: Main features of the Lk selected stacking bins: (1) bin number; (2) number of stacked sources; (3) average 
rest frame Lk in the bin (with corresponding standard deviation); (4) average redshift in the bin (with corresponding 
standard deviation); (5) average stellar age in the bin (with corresponding standard deviation); (6) average stellar mass 
in the bin (with corresponding standard deviation); (7) signal to noise ration of the stacked signal; (8) average 0.3-8 X- 
ray luminosity after LMXB subtraction estimated from a thermal model. The first uncertainty reported in parenthesis 
is due to standard deviation of counts, stellar age and stellar mass in the bin, while the second one takes into account 
the additional standard deviation of redshift.; (9) hardness ratio (with error due to standard deviation of counts, stellar 
age and stellar mass in the bin) 


(1) bin number 

(2) number of sources 

(3) < log(L^/L0) > 

(4) < z > 

(5) < log(age/yr) > 

(6) < log(M./Mo) > 

(7) SIN 

(8) < T-x-la/xb, 0.3-8 erg s ') > 

(9) HR 

1 

93 

9.38(0.22) 

0.16(0.04) 

9.66(0.26) 

8.93(0.27) 

3.12 

0.24(0.09,0.13) 

0.40(0.32) 

2 

946 

9.60(0.21) 

0.62(0.34) 

9.43(0.35) 

8.98(0.38) 

2.60 

11.40(5.07,10.72) 

0.24(0.38) 

3 

272 

10.60(0.22) 

0.27(0.07) 

9.60(0.24) 

10.07(0.33) 

3.09 

< 0.62 

> -0.45 

4 

190 

10.65(0.19) 

0.38(0.02) 

9.50(0.21) 

10.11(0.30) 

3.00 

2.11(1.45,0.26) 

0.37(0.61) 

5 

41 

10.67(0.19) 

0.44(0.01) 

9.43(0.19) 

10.11(0.32) 

3.01 

12.23(5.64,0.73) 

0.23(0.40) 

6 

43 

10.56(0.19) 

0.48(0.01) 

9.38(0.21) 

10.06(0.25) 

3.24 

21.48(8.62,0.67) 

0.61(0.33) 

7 

1017 

10.70(0.18) 

0.80(0.19) 

9.31(0.18) 

9.99(0.47) 

3.03 

34.53(19.32,8.27) 

0.40(0.39) 

8 

4 

11.32(0.15) 

0.18(0.01) 

10.01(0.03) 

11.01(0.17) 

3.14 

1.27(0.77,0.15) 

0.18(0.45) 

9 

7 

11.28(0.14) 

0.21(0.01) 

9.96(0.10) 

10.89(0.13) 

3.12 

1.57(0.87,0.06) 

0.05(0.52) 

10 

7 

11.17(0.15) 

0.22(0.01) 

9.86(0.23) 

10.78(0.15) 

3.02 

2.06(1.12,0.01) 

< 0.16 

11 

20 

11.43(0.17) 

0.24(0.02) 

9.94(0.09) 

11.05(0.20) 

3.24 

< 1.64 

< 0.02 

12 

43 

11.35(0.16) 

0.29(0.02) 

9.88(0.10) 

10.93(0.21) 

3.06 

1.07(0.96,0.13) 

-0.14(0.58) 

13 

59 

11.31(0.19) 

0.34(0.01) 

9.75(0.19) 

10.84(0.22) 

3.21 

1.71(1.60,0.13) 

< 0.15 

14 

13 

11.28(0.19) 

0.35(0.01) 

9.73(0.16) 

10.78(0.21) 

3.12 

7.73(3.44,0.18) 

> 0.06 

15 

13 

11.39(0.16) 

0.36(0.01) 

9.75(0.16) 

10.94(0.22) 

3.04 

6.53(3.61,0.01) 

0.16(0.47) 

16 

16 

11.43(0.21) 

0.37(0.01) 

9.72(0.22) 

10.91(0.25) 

3.04 

6.36(3.89,0.17) 

> -0.06 

17 

45 

11.37(0.18) 

0.38(0.01) 

9.68(0.22) 

10.88(0.23) 

3.05 

3.44(3.16,0.08) 

< 0.41 

18 

45 

11.55(0.22) 

0.39(0.01) 

9.70(0.19) 

11.02(0.28) 

3.03 

< 8.34 

< 0.31 

19 

91 

11.39(0.20) 

0.45(0.03) 

9.69(0.22) 

10.89(0.24) 

3.09 

6.92(5.57,1.03) 

0.23(0.76) 

20 

216 

11.44(0.21) 

0.59(0.06) 

9.61(0.25) 

10.92(0.28) 

3.04 

< 26.36 

-0.32(0.67) 

21 

27 

11.55(0.19) 

0.67(0.01) 

9.57(0.23) 

11.00(0.28) 

3.30 

45.92(19.33,0.08) 

0.36(0.45) 

22 

161 

11.46(0.21) 

0.69(0.01) 

9.56(0.23) 

10.91(0.29) 

3.06 

38.05(24.83,1.38) 

0.36(0.63) 

23 

238 

11.42(0.21) 

0.75(0.04) 

9.45(0.19) 

10.82(0.27) 

3.27 

45.23(30.76,5.62) 

< 0.25 

24 

59 

11.45(0.23) 

0.83(0.01) 

9.39(0.15) 

10.83(0.24) 

3.78 

111.59(38.68,0.45) 

-0.33(0.17) 

25 

319 

11.47(0.21) 

0.86(0.02) 

9.40(0.16) 

10.84(0.28) 

3.04 

65.09(47.91,3.10) 

> 0.11 

26 

122 

11.52(0.21) 

0.93(0.01) 

9.37(0.15) 

10.84(0.30) 

3.25 

104.73(50.43,2.76) 

-0.33(0.50) 

27 

565 

11.48(0.20) 

1.06(0.08) 

9.33(0.12) 

10.66(0.49) 

3.04 

204.40(159.37,38.54) 

< 0.14 

28 

24 

11.62(0.18) 

1.20(0.01) 

9.40(0.06) 

10.31(0.82) 

3.26 

1234.10(428.64,4.11) 

> 0.57 

29 

3 

11.48(0.12) 

1.21(0.01) 

9.40(0.03) 

10.85(0.14) 

3.09 

606.33(204.05,0.18) 

0.70(0.27) 

30 

112 

11.50(0.18) 

1.25(0.04) 

9.33(0.13) 

10.50(0.69) 

3.02 

653.77(277.84,46.88) 

> 0.31 

31 

80 

11.55(0.18) 

1.40(0.04) 

9.30(0.13) 

10.56(0.64) 

1.29 

< 1299.87 

> -0.89 

32 

8 

12.05(0.03) 

0.66(0.09) 

9.75(0.12) 

11.59(0.05) 

6.55 

368.43(62.19,123.45) 

-0.18(0.16) 

33 

27 

12.10(0.09) 

1.09(0.18) 

9.41(0.09) 

11.41(0.11) 

2.31 

332.98(211.58,137.06) 

-0.30(0.67) 






Table 2: Same as Table 1 but for redshift selected stacking bins. 


(1) bin number 

(2) number of sources 

(3) < log(L,r/L0) > 

(4)< z > 

(5) < log(age/yr) > 

(6) < log(M./Mo) > 

(7) S/N 

(8) < 7.x-LA/jfB,o.3-8/(10^ erg s ’) > 

(9) HR 

1 

186 

9.12(0.07) 

0,29(0.10) 

9.53(0.31) 

8.66(0.25) 

3.02 

0.96(0,36,0.72) 

0.54(0.33) 

2 

511 

9.84(0.23) 

0.31(0.10) 

9.52(0.26) 

9.28(0.34) 

3.03 

0.68(0.28,0.46) 

-0.33(0.52) 

3 

129 

10.42(0.06) 

0.34(0.10) 

9.49(0.24) 

9.96(0.12) 

3.07 

1.88(0.85,1.20) 

> 0.17 

4 

241 

10.75(0.10) 

0.35(0.08) 

9.56(0.22) 

10.29(0,15) 

3.01 

1.08(0.86,0.57) 

0.12(0.58) 

5 

97 

11.01(0.04) 

0.35(0.08) 

9.61(0.23) 

10.57(0.09) 

3.11 

1.77(1.33,0.94) 

> -0.34 

6 

9 

11.09(0.01) 

0.32(0.08) 

9.76(0.23) 

10,68(0.09) 

3.38 

6.38(2.36,3.77) 

-0.04(0.36) 

7 

91 

11.17(0.04) 

0.36(0.07) 

9.70(0.22) 

10.74(0.09) 

3.06 

2.38(1.84,1.03) 

< -0.12 

8 

37 

11.27(0.02) 

0.36(0.07) 

9.74(0.20) 

10.85(0.08) 

3.06 

4.25(2.42,2.08) 

0.07(0.48) 

9 

46 

11.35(0.03) 

0.35(0.07) 

9.81(0.19) 

10.95(0.08) 

3.03 

2.85(2,06,1.30) 

< -O.IO 

10 

25 

11.43(0.02) 

0.37(0.07) 

9.80(0.20) 

11.02(0.08) 

3.02 

6.17(3.61,2.70) 

> 0.42 

11 

20 

11.48(0.01) 

0.38(0.07) 

9.82(0.17) 

11.09(0.07) 

3.10 

8.66(4.38,3.50) 

-0.37(0.50) 

12 

18 

11.52(0.02) 

0.36(0.07) 

9.84(0.13) 

11,10(0.06) 

3.25 

5.70(3.02,2.77) 

-0.07(0.42) 

13 

28 

11.59(0.03) 

0.36(0.07) 

9.86(0.10) 

11.20(0.05) 

3.05 

3.96(2.80,1.82) 

< 0.50 

14 

13 

11.69(0.02) 

0.37(0.06) 

9.85(0.16) 

11.27(0.10) 

3.13 

8.90(5,04,3.12) 

< -0.33 

15 

18 

11.81(0.05) 

0.38(0.06) 

9.83(0.07) 

11.37(0.09) 

3.05 

6.00(4,60,2.10) 

-0.01(0.53) 

16 

425 

9.83(0.24) 

0.73(0.14) 

9,38(0.37) 

9.09(0.46) 

3.02 

43.08(15.98,20.05) 

0,52(0.36) 

17 

1159 

10.91(0.18) 

0.77(0.13) 

9,33(0.17) 

10.26(0,38) 

3.07 

21.41(15.95,9.18) 

0.05(0.54) 

18 

262 

11.28(0.04) 

0.79(0.13) 

9.41(0.20) 

10.74(0.15) 

3.03 

58.62(32.17,23.11) 

-0.33(0.42) 

19 

34 

11.36(0.01) 

0.82(0.14) 

9.45(0.21) 

10.79(0.21) 

3.30 

131.37(47.08,52.96) 

0.11(0.33) 

20 

111 

11.40(0.02) 

0.81(0.14) 

9.44(0.21) 

10.86(0,11) 

3.04 

103.10(48.81,43.24) 

0.20(0.42) 

21 

106 

11.47(0.02) 

0.81(0.14) 

9.45(0.21) 

10.93(0,14) 

3.08 

115.21(55.10,47.07) 

0,50(0.47) 

22 

274 

11.61(0.05) 

0.80(0.14) 

9.53(0.22) 

11.08(0.11) 

3.07 

52.12(37.39,21.70) 

> -0.29 

23 

8 

11.72(0.01) 

0.77(0.10) 

9.53(0.15) 

11.22(0.12) 

3.17 

204.51(79.18,65.47) 

0.18(0.37) 

24 

51 

11.75(0.02) 

0.79(0.15) 

9.60(0.24) 

11,24(0.10) 

3.78 

102.39(36.38,46.82) 

< -0.59 

25 

86 

11.84(0.04) 

0.81(0.14) 

9.59(0.20) 

11.32(0,10) 

3.03 

118.48(66.52,48.44) 

< 0.09 

26 

28 

11.95(0.02) 

0.85(0.12) 

9.57(0.19) 

11.41(0.09) 

3.24 

149.06(60.70,33.54) 

0.04(0.39) 

27 

15 

12.04(0.02) 

0.81(0.16) 

9.60(0.21) 

11.50(0.10) 

6.03 

397.82(77.91,190.80) 

-0.30(0.20) 

28 

8 

12.13(0.03) 

0.86(0.16) 

9.52(0.17) 

11.50(0.09) 

1.20 

< 269.42 

-0,16* 

29 

737 

11.14(0.27) 

1.18(0.12) 

9.30(0.20) 

9.93(0.79) 

3.06 

402.00(188.34,102.88) 

> 0.20 

30 

127 

11.68(0.06) 

1.20(0.12) 

9.36(0.10) 

10.90(0.44) 

3.25 

404.13(176.27,100.36) 

-0.23(0.44) 

31 

14 

11.84(0.02) 

1.25(0.15) 

9.34(0.11) 

11.19(0.16) 

3.18 

299.63(107.66,90.89) 

-0.46(0.41) 

32 

33 

12.01(0.13) 

1.21(0.13) 

9.38(0.05) 

11.24(0.41) 

2.32 

557.18(342.18,144.63) 

< 0.36 
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Fig. 3.— X-ray luminosity (after LMXBs subtraction) versus K-band luminosity for the stacking bins selected in Lk 
(Left Panel) and redshift (Right Panel). The horizontal and vertical black error bars indicate the statistical uncertainties 
on Lk and Lx-lmxb, respectively, while the vertical green error bars show the uncertainty on Lx-lmxb due to redshift 
spread in each bin. The stellar mass M* from COSMOS catalog is plotted on the top axis. The local relation by 
KF13 is represented with a full black line. The dashed lines are drawn to include the most Lx,gas luminous source 
ETG (M87) and the most Lk luminous ETG (NGC1316) in the Local Sample (BKEll, KE13), respectively. Bins are 
labelled as in Table 1 (for the left panel) and Table 2 (for the right panel): red. X-ray over-luminous with respect to 
the local Lx^gas - Lk relation; yellow. X-ray under-luminous bins; grey and green bins following the local relation with 
Lx-lmxb lower and higher than 10'** erg s *, respectively. Triangles pointing down represent upper limits. 
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Fig. 4.— Left; HRs of the L/f-first stacking bins versus redshift. Right: HRs of the z-first stacking bins versus 
Lf^. Triangles pointing up and down represent lower and upper limits on HR, respectively, while the box for bin 28 
represents unconstrained HR. 
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Fig. 5.— Same as Figure 3, but with bins labelled according to their average redshift, and HRs indicated by the slope 
of the orange segments. Black error bars reflect statistical uncertainties, vertical green bars show the uncertainty on 
Lx-lmxb due to redshift spread in each bin. Red boxes show the BKFl 1 and KF13 samples, blue circles show the C14 
sample of C-COSMOS X-ray detected ETGs. Triangles pointing down represent upper limits. 
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Fig. 6.— Left; HR-z scatter diagram of L/^-first bin. The vertical error bars indicate uncertainties on Lx-lmxb 
due to counts and to LMXB subtraction. Triangles pointing up and down represent lower and upper limits on HR, 
respectively. Squares represent unconstrained HRs. Dashed lines show HRs expected from thermal models (APEC) 
with increasing temperatures from bottom to top 0.7, 1, 2 and 3 keV. The full lines represent simulated HR from power- 
law models with hxed slope T = 2 and increasing intrinsic absorption form bottom to top 0, 10^*, 10^^ and 10^^ cm“^. 
Right; Same as left panel, but for models comprising a power-law PL with fixed slope F = 2 and intrinsic absorption of 
10^^ cm“^, plus a thermal model with temperature of 1 keV (APEC). The relative normalization of the two models is, 
from bottom to top, 10^^, 10^^, 10 *, 1, 10, 10^ and 10^, corresponding to a contribution of the power-law component 
to the total flux of ~ 0.2%, ~ 2%, ~ 15%, ~ 62%, ~ 94%, ~ 99% and 100%, respectively. 
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Fig. 7.— X-ray luminosity (after LMXBs subtraction) versus black hole mass Mgu (evaluated from M, according to 
the Kormendy & Ho (2013) relation) for the L/^-first stacking bins. The stellar mass M* from COSMOS catalog is 
plotted on the top axis. Lines of constant LboilLEddi given Mbh and Lt,oi, are shown as diagonals (from bottom to top 
10^^, 10^^, 10^^, 10^^ and 10^^). The bins are color-labeled as in Figure 3. Light blue boxes represent the predictions 
of the Volonteri et al. (2011) models for M, = 10® Mq, 10*°Mq and lO" Mq, where the top of each box refers to 
radiatively efficient accretion flows and the bottom refers to RlAFs. 












